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A direct calculation of the elements of the photon polarization vector for arbitrary momentum 
in the helicity basis shows that it is not a vector but a complex bivector. The bivector real and 
imaginary parts can be directly equated with electromagnetic field amplitudes and the associated 
field equations are the Maxwell equations in time-imaginary space. The bivector field exhibits 
a phase freedom (Berry, or geometric phase) dependent on the rotation history of the field or 
observer. Phase freedom is not intrinsically present in the longitudinal excitations of the field and a 
general argument connects quantization of angular momentum with the observation of phase changes 
associated with frame rotation. Current and translation operators can be defined for bivector fields 
that are free of defects associated with a quantized vector potential. 
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I. INTRODUCTION 



o 
o 

<N 

Wigner first pointed out that quantum states should be representations of the inhomogeneous Lorentz group 0|. 
However application of this simple intuitive principle has proven surprisingly challenging, particularly in the context 
of the massless spin one field 2}. One of the most widely encountered representations of the photon field is of a 
■ quantized Gibbs (column) vector potential 



A»(x) = (27T)- 3 / 2 [ 4=K(P, <r)o(p, oV px + e^*( P , ^(p, n)e^] (1) 

in which e A '(p, a), the photon polarization vector, satisfies p M e' i (p, a) — 0, e°(p, a) — 0. For a pure motion in the 
1-direction 

e(p,±(r) = (0,0,±* > l)/V5 (2) 



o 
o 

Because the polarization vectors are transverse, they impart to the field the influence of holonomy, a quantity associated 
with the change in orientation of a tangent vector directed by parallel transport around a non-Euclidean surface. A 
substantial literature has developed around the discovery by Berry 3] that the transport of a quantum system through 
an adiabatic circuit induces a change in phase that is dependent on the area of the surface enclosed and the angular 
momentum of the system. The identification of this effect with holonomy Q , its realization experimentally through 
experiments with light 0, [|| , the elucidation of its origins in the Poincare group @, II, El , and the ensuing debate 
over whether the experimentally observed effects have classical or quantum origins |lCt 111! ITJ . Il3| have recently 
been reviewed comprehensively jlif . Berry's original treatment has been extended to non adiabatic circuits 0] and 
open trajectories 0,EJ; and the contributions of holonomy were recognized to be implicit in earlier studies 
of polarization phase by Pancharatnam [l9l |20) and in the proposal by Aharanov and Bohm |2l| to measure phase 
changes in non-simply connected topologies. 

In this work we carry out a direct calculation to demonstrate that the photon polarization 'vector' is not a vector but 
a complex bivector, and that the bivector can be equated with the electromagnetic field, not the vector potential. The 
bivector field exhibits a phase freedom associated with rotation history that cannot be separated from the observer 
frame. If the widely used Gibbs vector formalism is extended to spaces in which the time dimension has imaginary 
weight, a covariant Hamiltonian can be created that yields well-behaved eigensolutions with positive and negative 
energy eigenvalues and that yields the Maxwell equations upon application to the fields. Current and translation 
operators can be defined that are free of various defects associated with the quantized vector potential. Experimental 
tests can be devised that distinguish between a vector potential and a bivector, using a simple idealization of a single 
mode state to explore contrasting predictions of the dependence of the energy of a photon on its modulation envelope. 



II. A COVARIANT OPERATOR ALGEBRA FOR GIBBS VECTORS 



In the following the metric signature is ( — h H — h) and 'a pure motion in the n-direction' signifies that the velocities 
along the spatial axes j ^ n are null. Our first task is to develop a covariant Gibbs vector algebra that allows a 
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Hamiltonian equation for the electromagnetic field. To achieve this we introduce operators x M thai can be considered 
analogs, in a vector representation, of the Weyl operators in the spinor representation. The components 
(formally, X^p) are rank two tensors that can be written 
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(3) 



The x M % an d the related X^v 9 constitute two subalgebras of Mai(4, C) that are isomorphic to the left and right 
subalgebras of the real matrix representation of quaternion multiplication Mat(4,M.) ~ Hcg)H The three spatial 
operators x 3 have eigenvalues of ±1 and eigensolutions that are irreducible representations of the Lorentz group 
in a space that assigns imaginary weight to the (scalar) time direction. In this space the boost transformations 
are reflections, so the fundamental group operation is a reflection. Reflections are in any case more primitive than 
rotations since by a theorem of Hamilton any rotation can be decomposed as a pair of reflections, e.g. .23]. The 
normalized eigensolutions of xi are: 




(4) 



with respective eigenvalues {1,-1,1,-1}. The first pair are eigenbivectors or simply bivectors and correspond to the 
two orthogonal states of the polarization 'vector' JSJ; the second pair are time-imaginary vectors. Like the a 3 , the 
Xj , j S 1,2,3 form the basis of a Clifford or geometric algebra CI3 on M 3 defined in the usual way by 



XiXj +XjXi = 2$ijU 



(5) 



where I4 is the 4x4 identity matrix. By J5J) appropriately normed Clifford numbers composed from the x % are 
elements of the (odd + even) group Pin(3) with Pin(3)/{±1} ~ SO(3), since Pin(n) forms a double cover of SO(n), 
[22ll24] |. The x % can be divided into spatial and temporal parts x % = K l + S l = K i + S where 



Si u, — Si ^0 — 



Si 



(6) 



K o c 
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(7) 



The Si are the generators of the spin one rotation operator about axis i, often written ^fc, a nd have, with the Ki, 
the following properties: 



S l 2 + Ki 2 



Sj h S i k l +K i 1 k K i k l = 6S l 



(8) 



e KS * = K 2 + Si 2 cosh k + Si sinh re (9) 

The Si are also the generators of bivector boosts in direction i. The Ki are the generators of vector boosts in direction 
% and satisfy 

e K K * = S 2 + K 2 cosh K + Ki sinh re (10) 

The first two and last two elements of Q are eigen(bi)vectors of Si and K\ respectively. 
The Si and Ki obey the commutation relations 

[Si 1 Sj ] — i tij fc Sk 
[Ki,Kj] ^itijkS k 
[S i} Kj] =it tjk K k 

and the above relations differ from the usual formulation only in the sign of the righthand side of the middle equation. 
However as a result the operators Xi conform to the Lie algebra of the quaternions, 

[Xi,Xj] = ZieijkXk (11) 
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The Hamiltonian operator x J Pj has doubly degenerate eigenvalues ±po and hence admits a covariant equation of 
the form 

X%V = (12) 
as well as a contravariant (opposite helicity) equation: 

XW = (13) 

In the case of a pure 1-direction motion the first eigenbivector of (@J, V+ — (0, 0, — i, l)/v2 is a solution of l|12() 
and the second, ^_ = (0,0, i, l)/\/2 is a solution of (jT3|l . and these correspond to the positive energy solutions for 
a positive and negative helicity particle, respectively. To the ip± can be associated complex conjugates i/j±; their 
modulus, or length, tp± ■ ip± is unity. Like their counterparts in the massless spin 1/2 held, tp + and ip_ undergo 
boost transformations by different laws, but behave the same way under rotations. The ip± can be rotated to 
an arbitrary direction of propagation by i\)± = e~ lX2 S2 e~ tXs S3r t/j± where A2 = arcsin(— p 3 j \Jp\ + p 2 ) and A3 = 
arcsin (p 2 / yV\ + P 2 + P% ) > the order of multiplication of the rotation operators being given by the leftmost first, 
rightmost last prescription (or in a spinor context, outermost first, innermost last), see e.g. |23l| sec. 39. This gives 

z 2 (p) := e -^e-^V + = ( °'^ lP2 ~ P0P3 ;" ?(Pl2+P32) ' P0Pl+?:P2P3) (14) 

v2po Vpi 2 +P3 2 

Similar actions on pure 2- and 3-direction bivectors, <ft+ — (0, 1, 0, — i) /a/2 and = (0, — i, 1, 0)/\/2 respectively, 
lead to 

»(p) := e-^e-^ + = ( °' P ° P2 + ' PlP3 ' ~ P ° Pl + ™' (V + ^ (15) 

where ^3 = arcsin(— p\j ^Jp\ + p\ and fix = a,rcsin(p3 / ^/ p 2 +j»2 + P3) and 

z M ■= e -*** e -*>*£ = ( °' ~ l ^ 2 + P3 ^ ' lpiP2+PoP3 > - po + » Pi P3) 

V2P0 V(P2 2 +P3 2 ) 

where v x = arcsin(— p 2 / \/p\ + p 2 ) and t/ 2 = arcsin(pi/y / pf _ T^Tpf). A structure closely related to z 3 can be found 
as a column element of a unitary matrix governing transformation of spin one states Q • 

The initial values of ip+ , and £ + above have been chosen for consistency under rotation of the coordinate axes 
between the positive helicity bivector solutions and the positive momentum vector solutions 

(i ) l,0,0),(i > 0,l,0),(i, 0,0,1) (17) 

The forward cycle (1 — > 2 — > 3 — > 1) of transformations of both the bivectors and the vectors are produced by 

successive action of the axis rotation operator e 3 ^ . The reverse cycle (1< — 2< — 3< — l)is generated by 
the operator formed from the inverse (= the square) of the positive cycle operator. 

Despite their apparent dissimilarities, each of the representations i|14l - I16fl is equivalent; each has unit length, 
gives rise to energy-momentum currents of the same magnitude (discussed below) and transforms under boosts and 
rotations in the same way. However each Zj is poorly defined for a pure motion in direction j. The normalized sum 
of bivectors 



Vpz 2 + P3 2 zi(p) + \Jp\ 2 +P3 2 z 2 (p) + \/pi 2 +P2 2 z 3 (p) 

zs{p) := / 2 2 2 ' ' 

V {Pi ~ P2) 2 + (P2 - P3) 2 + (-Pi + P3) 2 



with components 



zs(p)o = 

£fc = i {-pgtjkm p m +i (pj Pk-Pk 2 )) 
V2po \/ (JP1-P2) 2 + (pi-ps) 2 + (P2-Ps) 



Zs (pj . _ Lk=i{-Po £ jk m p m +t(Pj Pk-Pk )) ^ ig ^ 



(with implicit summation over m) is free of this defect. The subscript S reflects the space-like character of the 
polarization vector. 
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Upon rotation to arbitrary momentum the negative eigenvalue solutions in yield the complex conjugates of i|14l 
- I16|) and (|18p. corresponding to the solutions for negative helicity and denoted by overlining. Because 

z s z s = z s z s = 1 z s z s = z s z s = (20) 

the 2s and zs recapitulate the familiar orthogonality relations of the polarization vector for distinct helicities a, a 1 : 

e%(p, ( r')e' I (p,a) =5^ (21) 

In a similar manner we find that 

^T f e li (p,CT)e* v (p,a)= z s ®z s + z s ® z s 

cr 

= otherwise (22) 

which is the standard relation for the projection matrix of the polarization vector. 



III. IDENTITY WITH ELECTROMAGNETIC FIELD 

The Zj are solutions of i|12|) . but describe a field that is not a simple vector potential. Using the correspondence 
between the symmetric field stress-energy tensor O Q/3 := —g QAI F x/3 — \ g Q/3 F /iA and the stress-energy tensor 
T Q/3 with components p a pp/po to translate between momentum and field representations, we find 

Zj { P ) = _ZL ~*h^J L(0,ibi + e 1 ,ib 2 + e 2 ,ib 3 + e 3 ) (23) 

and 

z s (p) = =J= lbj+€j (0,ib 1 + e 1 ,ib 2 + e 2 ,ib 3 + e ;i ) (24) 

^ A /(E,^) 2 + (E, 



The differences in representation reflect, as in the momentum formulation, alternative normalizations differing in 
phase. The normalization factor of l/y/2po emerges naturally as a consequence of the unit length of the bivector. 

With H23|) and l|24() we have found that explicit application of the program to represent particles as irreducible 
representations of the inhomogeneous Lorentz group in the Gibbs vector representation leads to the conclusion that 
the elementary excitation of the massless spin one field is not a vector but a complex bivector. In the normalized 
form of (|18|l . the real and imaginary parts can be equated with the electric and magnetic amplitudes of the photon 
field respectively. That is, (|18|1 can be decomposed as 



where, for example, 



zs(p) ■= ze(p) + i zb(p) (25) 

ze(p)i = po(p 2 -ps)/N 

zb(p)i= (-Po 2 +PiEi=ift)/^ 

N= V2po sj (pi - p 2 f + (pi - Psf + (p 2 - Pa) 2 (26) 

ze(p) reverses upon spatial reflection (pi — > —pi) but not time reflection (pa — > — po), whereas zb(p) shows the opposite 
behavior, consistent with the known symmetries of the electric and magnetic fields. The ze(p) and zb(p) constitute 
an orthogonal basis from which the projection matrix can also be derived. 

zs ® zs + zs ® zs = 2 (ze(p) ® ze(p) + zb(p) ® zb(p)) (27) 
With the results above we can define a positive helicity spin one field 

ip+ = —=={Q,ibi + ei,ib 2 + e 2 ,ib 3 + e 3 ) (28) 
V2p 
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in which the factor e 10 assimilates the phase freedom of the field. The field equations 

e 10 e iB 

X M d^+ = (-V-b + iV.e,a e-Vxb + !(9 b + Vx e ))^= = (*P,-j)-«= ( 29 ) 

y/2p V 2 Po 

are related to the classical Maxwell equations in a time-imaginary space. The negative helicity field, which is the 
complex conjugate, follows the equation 

X »d^ = (tp,j)-^= (30) 

so the conventional representation, is not a solution of a single field equation and the field must be represented 
instead as a direct sum ip+ © . Application of the positive helicity equation to the negative helicity field or vice 
versa results in a mispairing of the sign of the time derivative with that of the cross product, i.e. yields equations of 
the form 

<9e + V x b = and <9b-Vxe = (31) 
The sourceless field equation associated with the direct sum can be written 

*(£ 4^) (-&)-* (40-° < 32 > 

Qq being the operator inner product, an 8 x 8 matrix. 



IV. TRANSFORMATION PROPERTIES OF THE FIELDS 



The boost transformations of zs and zg differ 

z s (p') = e- KS *z s (p) = z s (p)e KS > z s (p') = e« s *z s (p) = ^(p)e" KS< (33) 
whereas the rotations are the same 

z s (p') = e iK3 *z s (p) = z s (p)e-^ s > z S {p') = e lKS ^z s (p) = zsW^ 3 * (34) 
Applied to the representation l|32|l this implies that the generators of boosts and rotations have the forms 

(£-*) ™ d (£S) (35) 

respectively. Although the Gibbs vector formalism is the most frequently encountered depiction of the electromagnetic 
field, the most compact way to represent the field operators is not by Gibbs bivectors but by spinor bivectors (or, 
more directly, by elements of a Clifford algebra |22l l25l l2a |). which will be presented elsewhere. The time and space 
reflection operators defined by 

*=-0u"fc) ™ d Ji) < 36 » 

have action on the operator Qq of (|32|l given by the inner products RtGoRt and RsQdRs respectively, consistent 
with their character as reflections, not rotations. 



V. PHASE FREEDOM AND TOPOLOGY 



Because the lengths of the bivectors (|14l - I16H are equal, the differences between them constitute variations in 
normalization that affect the argument but not modulus of a complex quantity, a form of phase freedom that can be 
identified with the well-studied Berry phase 0, 0, 0] , recently reviewed in [ljj . The phase has a history dependence 
which is easily displayed without reference to circuit as follows. 

If the order of rotations is reversed for each of the (|14l - I16fl . a different collection of bivectors results. For 
example, from the pure 1-direction bivector %p± we can carry out the rotation ip±' = e^ 1 ^ 3 Ss e~' lX2 S2 ip± where 
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A2 = arcsin(— p 3 /\/Pi + P% + P3) an d ^3 = arcsin(p 2 /\/Pi +p|)- ^ n ^ s case the resulting is poorly denned 
for motion in the 3-direction, not the 2-direction as in 1141) . The collection of bivectors under reversal of the rotation 
order then are: 



0+ -> 2/3, C+ S/i, f+ -> y2 (37) 



with the translation to field components 



Vj = -i= fi^ L(0,ibi + e 1 ,ib 2 + e 2 ,ib 3 + e 3 ) (38) 



differing from l|23|) by a factor of i. Thus we can easily calculate changes in phase without completing a circuit. 
Alternatively, any open path can be turned into at least one closed circuit by joining the endpoints by a geodesic, 
along which the angle of rotation by parallel transport is zero 0, 0] 

By rotating from a pure motion in the 1-direction to a motion in an arbitrary direction via the sequence leading to 
(|23|l , and then back to a pure motion in the 1-direction by the reverse of the sequence leading to (|38|1 , a rectilinear 
path is inscribed on a sphere, and the regenerated bivector at the close of the loop has the form 

zbt arctan( P2P3 ) , /on\ 
e y POPl'ip ± (39) 

The simplest interpretation of 1)39(1 is that the positive helicity bivector has been rotated by arctan(^2.) and the 
negative helicity bivector has been rotated by — arctan(|^-). However if a positive helicity state is rotated about the 
axis of propagation by an amount 7 and a negative helicity state is rotated about the same axis by —7, the effect of 
rotation is equivalent to a translation in time of the phase of both states by the same amount. In particular, for any 
equally weighted mixture of such states the holonomy angle should be zero. This conclusion would be consistent with 
the inferred dependence of phase on angular momentum jj, 0, , but leads to the prediction that a linearly polarized 
beam, which has obvious transverse character, should show no holonomy. This is not observed experimentally [||. 
Although it might be possible to invoke the influence of optical fiber geometry on rotation of the polarization plane, 
considerations based on the orientation of the field with respect to time allow an alternate interpretation. 
The classical sum of helicity states described in the preceding paragraph is |26j chap. 7 

e iv^"\ + ) + e -ip^\-) (40) 

where for avoidance of doubt we can multiply both states by the same modulation envelope, such as a simple Gaussian 
functional (2na 2 /Po) -1 / 4 exp(— (£V a;AI ) 2 /(4c 2 )) H51- and confirm that the resulting structure describes a plane wavelet 
propagating forward in time. (|40ll has no obvious pathologies from a classical perspective. However it docs not allow 
a Hamiltonian equation that is consistent with all properties of the field, even if the sum is taken to be a direct tensor 
addition, because the derivative with respect to time yields both positive and negative energies and the resulting 
currents are then pseudovectors, not vectors. Consistent with this, if l|4U|l is multiplied by a Gaussian modulation 
function the Fourier transform of the resulting state is a sum of two momentum amplitude components, one centered 
on p and the other on —p. If we require a common orientation for all states in both time and time derivative, we have 
a direct sum of the form (|32H . and the Hamiltonian has the form 

*-(££)-C8F-£,) <41) 

which has quadruply degenerate eigenvalues of ±Po- The time evolution of states is given by 

' e TW) (-£.)«*- =(-1) e,p "" <«> 

With this, multiplication of ip+ and ip_ by e %1 and e~ lJ , respectively, is equivalent to phase translation in opposite 
time directions, leading to rotation of the plane of polarization. Representing the polarization state by the parameters 
of the Poincare spinor (helicity basis), we have for the positive energy part of the field 



i>+ a+ 



e ip » x " (43) 



where a+ and a_ are the annihilation operators for the positive and negative helicity states, respectively. In the 
standard depiction, the 'north' and 'south' poles of the Poincare sphere are pure + and — helicity states represented 
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by amplitude weights of cos 0/2 and e~ z ^ sin 0/2 respectively, and the equatorial antipodes (front and back) are pure 
'x' and 'y' states. The plane 8 = ir/2 is the plane of linear polarization, and 8 = 7r/2, (j> = corresponds to linear 
x axis polarization while 8 = n/2,(f) = n corresponds to linear y axis polarization. From this we can see that the 
effect of multiplication of an equally weighted mixture of positive and negative states (8 = ir/2) by e 11 and e~ 27 , 
respectively, is equivalent to a rotation of the plane of linear polarization by 7; a corollary is that the phase factors 
of (Kfflj) do not have a simple geometric interpretation. With respect to the angular momentum dependence found by 
Berry, we must also conclude that the equally weighted direct sum of time-oriented (1,0) and (0,1) representations 
does not transform as a state of zero angular momentum under rotations. 

In the standard spherical coordinates p\ — po cos(/>sin#,p2 = Po sin0sin#,p3 = po cos 8, (|39|l gives 

g=f=* arctan(cos 9 tan <f>) (44) 

and by the Gauss-Bonnet theorem, for a large class of simple loops that enclose an area that can be represented as 
the union of a countable collection of triangles, the exponent is equal to the area of the spherical surface enclosed 
divided by the square of the radius of curvature, \p\. For the calculation above this is easily confirmed in the limit 
8 — > 0. The general result can be reached in a direct and elegant way from the properties of the holonomy group and 
the unitary representation of the massless spin one little group 0. 

Reversal of the direction by which the path is traversed leads to reversal of the sign of the exponent. In addition it 
can easily be seen that for a spin n field, constructed as the outer product of n spin 1 fields, the phase acquired upon 
traversal of the circuit would be n times the value obtained upon traversal by a spin one field (apply the sequence of 
transformations to each element of the outer product). The related SU(2) phase value ±arctan(cos6?tan^) [27l |2S| 
has been identified with the Pancharatnam phase |19| , obtained by rotating the polarization state about a closed path 
on the Poincare sphere. 

Rotation of the vector solutions of (QJ to an arbitrary direction of motion by recapitulation of the steps leading to 
(|14l - ll6fl is simple and yields the same result in all three cases, 

z t{p) ■= -7= — (ip ,pi,p 2 ,P3) (45) 
V2po 

the longitudinal field, with time-like polarization. There is no phase freedom in this case because the solutions 
describe polar, as opposed to tangent, vectors; there is correspondingly no quantization known to be associated with 
the longitudinal field. 

The prescription for boost transformation of zt differs from that for z$, and as was the case for the zs, the 
transformation laws for zt and ~zt are different: 

z T (p') = e KK >z T (p) = z T (p)e~ KKi z t (p') = ^ kK ^t{p) = z T (p)e K ^ (46) 

The rotation operator is the same for zs and zt (and "zs and zt) 

z T (p') = e iK ^z T (p) = z T ( P )e~ iKS i (47) 



VI. PHASE FREEDOM AND QUANTIZATION 



Phase freedom imposes additional constraints on field structure. Because the wave-like behavior of elementary 
particles is experimentally well-validated there must be some principle that prevents the phase of the field from 
changing arbitrarily from one space-time coordinate to another. If this were not present all phases could equally 
coexist at all space-time coordinates - or equivalently be unpredictable - and interference phenomena would not be 
observable. The simplest reconciliation of phase freedom with wave-like behavior is to posit that the phase of a photon 
is determinate, but changes upon rotation from one coordinate system to another in a path-dependent manner. That 
is, that phase is embedded topologically and depends on the rotation history of both the particle and the observer. 

Now imagine that the photon field were not quantized, specifically that the angular momentum could be varied 
continuously by some unspecified classical mechanism. The angular momentum in this hypothetical field would 
presumably be represented by an intensive variable, i.e., would be localizable. Since the phase change associated with 
variation in path orientation depends linearly on angular momentum, frame rotation of a wave train would result 
in a phase change proportional to the total angular momentum accumulated at the point of measurement, in turn 
dependent on the length of the electromagnetic impulse measured. We could formalize this by saying the field would 
have an additional multiplicative factor e lx ° e where 9 is proportional to the solid angle of the circuit traversed (in the 
case of a closed loop). However this would be equivalent to a change in energy of the system by — HO. If this were 
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not sufficiently disquieting, by the analysis above, the plane of linear polarization would rotate with time. Although 
it has been argued [TTL Il2| that the rotation of the plane of polarization of linearly polarized light injected into a 
helically wound optical fiber represents a classical effect, the necessity for a discrete angular momentum, and 

hence field quantization, appears unavoidable. 



VII. CURRENT AND GRADIENT OPERATORS 



For massive and massless spin 1/2 fields, there exist well behaved four-vector current operators whose expectations, 
typically written tp^^ip and ip^a^ip, respectively, yield the predicted particle flux. If these did not exist there would be 
difficulties accounting for conservation of charge. Because photons do not bear charge, there is no obvious necessity for 
a conserved current, but there are related requirements to preserve conservation of energy- momentum. The standard 
definition of the quantized quantity corresponding to the Poynting energy-momentum flux vector is an expectation 
over a field and the time derivative of its adjoint: 

= - I d 3 xA^\x,t)d/dx^A^ +) (x,t) (48) 



This has been formally justified in the context of the canonical quantization of boson fields, which are thought to 
have canonical adjoints that are time derivatives of the ordinary adjoints; the canonical adjoint of spinor fields is the 
ordinary adjoint. Depending on one's perspective, either quantization, boson or fermion, could appear unnatural. For 
example, in the fermion case the canonical momentum is the adjoint field. 

From field equations of the form d fl x fJ " l P+ — Oj it would be reasonable to attempt to identify the x* with current 
operators, yielding operator expectations of the form Jd 3 i[f/' + 'x''^+] = j^, while allowing translation operator 
expectations to be formed in the usual way: J d 3 x [ip+^d^ V+] — P M - Although such a representation can be realized, 
the current operators transform not as vectors but as tensors, because the fields themselves are tensors. 

Our first task is to introduce an adjoint field as a formal conjugate, written in terms of the electric and magnetic 
amplitudes as 

V;+ f = -ts=(0, -i &i f + e x \ -i 6 2 + + e 2 \ -i 6 3 f + e 3 f ) (49) 

In the presence of interactions the adjoint field is an independent quantity and by convention labels the exiting particle. 
The momentum function conjugates corresponding to this field satisfy the adjoint equation z(p)p fJ ' = as well as 
the contravariant (negative helicity) equation p^x^^ip) = 0. The latter correspond to the Maxwell equations 

-i9 -id 

a^x^+ f = (V • b f -HV • e*, -doe* iVxb'l %{d^ + V x ^))—= = Up* J)—— (50) 

y/2p V 2 Po 

where we have superscripted the exiting current to signify its possibly different structure from that associated with 
the entering field. 

The Poynting currents for the field can be calculated as elements of the symmetrical tensor formed from the field 
and adjoint, expressed in components as 

^± tx x\y^±, (51) 

The currents formed from the amplitudes can be written 

V^VxVi = V'iVV'i = ±e ijk (e/ b k - 6/ e k qp i (&/ b k + e/ e fe )) /(2p ) (52) 

and 

ikfcVxVi = i^bj + e/e, ± i fa% ~ &/ ej ))/(2p ) (53) 

The spatial elements of the stress-energy tensor are 

^± < X l X m ^± = $lm{bj*bj + ejUj ± i (e/6j - b^e,)) 
-(h'b m + 6 m t 6/ + ejte m + e m f e/ ± i(e^b m - 6 m »ej + e m f 6/ - 6/ f e m )) (54) 
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In the absence of interaction the field component operators b%* , ej etc., must form Hermitian products (i.e., bi'ej = 
ej'bi) to insure that the momenta, etc., are real. With this the imaginary terms in (|52I54|) vanish and 

iktVVi = ±e<jfc (e/ b k - e fc t fy) /(2p ) (55) 

and 

^± f X V± = + e/ ej )/(2p ) (56) 

which reduce to the classical expressions. The Hermiticity constraint does not apply in the presence of interactions, 
as the incoming and outgoing fields have different energy-momentum structures. 

In the absence of interaction the current bilinears generated from the momentum functions of the fields and adjoints 
have the properties of a generalized velocity, i.e. 

zt{p)X^t{p) =zs{P)X^s{p) =P»/p° (57) 

although it should be kept in mind that in the case of the zs, the currents are components of a tensor, not a vector. 
For example the full relation for the stress-energy tensor is 

zs(p)x a x^ v zs(p)=P a pP/pl (58) 

and the action of the boost operator on the currents is a tensor reflection, 

jn> = ^e KS X fl e KS ip (59) 

whereas the action of the rotation operator is a tensor rotation 

j»> =^ e - iKS x^e iKS ip (60) 

The boost operation is a reflection because the current operator corresponding to the time dimension is a scalar; if it 
were a vector element (the usual case) the boost operation would be a rotation. 

The field expectation of translation operators is, by most definitions, computed over the scalar part of the field. In 
the case of a scalar particle, there is no current operator known apart from the translation operator, and the definition 



p" = J d 3 x : [4> + % i>+ - ^+d, : (61) 

where the colons denote normal ordering of creation and annihilation operators, can be identified with the operator 
expectation. This definition extends to both spinor and bivector fields, and to higher order fields for which the adjoint 
can be identified as the field quantity that yields a probability density upon contraction with the field. 



VIII. ORTHOGONALITY RELATIONS 



The zt and zs are orthonormal with respect to inner product (with the application of the energy-momentum 
constraint p^p^ — 0): 

zt zt — zs z s = 1 

z T zs = z~s zt = (62) 

and the electric and magnetic amplitudes in the transverse representation have the additional normalization and 
orthogonality properties (with the energy-momentum constraint) 

ze(p) zb(p) = 
z T (p) ze(p) = z T (p) zb(p) = 

ze(p) ze{p) = zb{p) zb(p) = 1/2 (63) 

consistent with the idea that an orthogonal basis can be formed from vectors parallel to the momentum, the electrical 
field and the magnetic field, and that the electric and magnetic components of a source-free simple field contribute 
equally to the energy density of the field. 
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IX. PROJECTION OPERATORS 



Projection of positive and negative energy states requires an idempotent operator, usually formed as a sum over 
(orthogonal) spin states of the normalized outer product of basis vectors and their adjoints, P = ^2 a ip ® ip'. In the 
representation developed here the sum is not over spin states but rather longitudinal and transverse components. 

z.s <8> z s = Si j /2 - pi pj /2pa 2 - ie,; j k Pk /2po 1 < i, j < 3 

= otherwise (64) 

The spatial components of zt <8> ~zt are p% pj /2po 2 , leading to 

P + = z T <8> zt + z s <8> zs — — (65) 

2po 



P-=z T ®z T + z s ®z s = \ (66) 

2po 

These have the required properties P±P± = P± and P±P^ = 0. Also P+ z — z;P-~z = ~z; P + ~z = P_ z — 0. 
Because the negative energy solutions are indistinguishable from lefthanded solutions in the case of the transverse 
field, these operators also project out spin states. The structure of the projection operators illustrates the principle 
that both transverse and longitudinal contributions must be included to form a complete description of the field. That 
both contributions are required is also implicit in the description of the Maxwell equations. 



X. CONNECTION WITH OTHER FORMALISMS 



In some treatments the relation —dA/dt — E (in radiation gauge) is taken to justify the representation of the 
photon in terms of a quantized electric field 

/oo j 
d 3 pJ^-[e f "{p,a)a{p,a)e lpx +e' 1 *{p,a)a\p,a)e- vx } (67) 
-oo V 2 

in which the polarization vectors e M (p, a) are the same as those used to describe the vector potential. The magnetic 
amplitudes are given by 

H»(x) = (2n)-^Y r ^.fe (pXe(p ' ff))M fl (p,^ + (pXe(p ^V (p,.)e-n (68) 
a J-oo V 2 po po 

This depiction inherits some of the difficulties of the vector potential representation. The electromagnetic field 
transforms not as a vector Ij46fl but as a tensor the polarization bivectors e(p, a) describe the complex sum of 
electric and magnetic fields, not the individual fields in isolation; and the field conjugates transform by different laws 
and obey distinct field equations. 

The approach taken in this work appears to be most closely related to that of Bialynicki-Birula and Bialynicka- 
Birula [7j who have emphasized the role of the Poincare group in determining the structure of the fields. A gauge 
dependent covariant derivative operator does not feature in this analysis, however - the comparable action is fulfilled 
by operator inner products of the form I|12ll3fl in the momentum representation. 



XI. MODULATION ENVELOPE TERMS 



A vector potential with current given by JJSJl, and a bivector field with current based on l|52|) or l|61|l can be 
distinguished by the presence or absence of terms dependent on the modulation envelope. Introducing the generic 
real envelope functional /(p^a; M ), subject to the constraint J f 2 d 3 x — 1 we note that, as in the classical case [26^ . 
V* = f(p^x^)e WfiX,i is a solution of the wave equation d^d^tp = 0. Taking a pure polarization state, for example 
pure positive helicity light propagating in the 1-direction in a discrete single mode with wavenumber k, the vector 
potential wavelet 

A(x") = ty^=±{ e {k, <T)e lfc ^V + e(k, a)* e - lk » x "a^ := A + (x») + A~(x^) (69) 
V2wfc 
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gives, by (|48|l . a field momentum 

k A " = - [ d 3 x[d Q A-(x")d»A + (xn} =a k f a k / d 3 x[fc"(/ 2 + /'*)] = a k t ak + / d 3 x[f' 2 }) (70) 



whereas the bivector field wavelet 

B{x") = f{k^x^)(e(k, a)e lk ^"a k © ef^aj'e-*" 1 "^ := S+(x^) © B~(x^) (71) 
gives, by (|61|l . a field momentum 



J d 3 x: [B~{x v )d^B+{x u )~ B+{x u )d^B-{x v )] : = a^a k k'" (72) 



The difference between l|72|) and (f TCIf) is the spatial integral of a squared functional derivative. To see what this 
contributes we evaluate it for the case of the normalized plane wavelet envelope 

/fax") = (2^ 2 /fc 2 )- 1 /4 e -^j|^ (73) 

specialized to a pure 1-direction motion (integration over the 2- and 3- directions is taken to give unity) . In this case 
the Fourier transform in the two spatial variables x$, x\ of the positive frequency part of the field gives 

/(fepx^je**"*" e - ipoXo+plX1 dx d Xl = (87rcr 2 /fco) 1/4 e ^ ~ *(po-J>i) (74) 
and the squared derivative term in (|70|l gives 

A \f'h f A M (k XQ ~ k lXl ) 2 _ (fcQttQ-fcl*l) 2 1 

]=J dx x [ 25/2(j5nl/2 e ^ ] = ^ (75) 



so that the current (|70|l is 



k A ' 1 = (i k W(l + ^) (76) 



An energy-momentum contribution that increases with decreasing width of the modulation envelope should be demon- 
strable experimentally. 



XII. CURRENTS FROM FIELD SUPERPOSITION 



Another setting that allows discrimination of vector and bivector fields is the prediction of currents arising from field 
superposition. Consider for instance a field that is constituted by the division and redirection of two components. Let 
the two components be directed toward each other symmetrically after some experimental manipulation to separate 
them, and for convenience, let them be in the same polarization state. If the envelope functional of one component 
is f(k^x^) and of the other is g(-k^x^) (i.e., they are directed with equal and opposite momenta, evolving in the 
same time direction), and if we set the time of intersection of the fields to be the origin, the positive frequency part 
of the bivector superimposed field will be 

B + {x») = f{k^x»)e(k,<j)e lk » x " + gi-k^x^k.ay- 1 ^*" a* (77) 

in which the fields share the same creation and annihilation operators because they are derived from a single source. 
Using H77f) and (|61|l we find for the spatial components of the momenta 

<k B J >=a k t a k fc J J d 3 x[f - g 2 + 2Sin(2p ■ x)(f ' g + g'f)} (78) 

where the subtraction of squared functionals reflects their opposite contributions to total momentum; for the energy 
we find 



< k B ° >=a k t a k fc° / d 3 x[f 2 +g 2 + 2fgCos(2p-x)+2Sin(2p-x)(f'g-g'f)} 



(79) 
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so the energy is given by the sum of the energies of the components, as expected. For the case of a pure 1-direction 
motion, with a Gaussian envelope of the form l|73|) for /, and with g given by the same formula with opposite sign of 
pi, the interference term in (|78|) is given by 



2 2 i 2 2 

2 p x o sin(2 pi xi ) (80) 



V2 7TCT 2 CT 2 

Analysis of the currents by the vector potential treatment gives 

< k A " > = a k ta k fc" /" d 3 x[/ 2 + f' 2 - g 2 - g' 2 + 2iStn(2p ■ x)(/'.g + f g') + 2 i Cos(2p ■ x)(/.g' - f'g)} (81] 

hence in addition to the squared functional derivatives encountered in l|7()|l there are imaginary interference terms. 
I thank Jeremy Gunawardena, Tom McCauley, and Rebecca Ward for suggestions and encouragement of this work 
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